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Abstract 

A new relaxation mechanism is shown to arise from overdamped two-level 
systems above a critical temperature T* ~ 5 K, thus yielding an explanation 
for experimental observations in dielectric glasses in the temperature range 
between T* and the relaxation peak at 50 K. Using the distribution function 
of the tunnelling model for the parameters of the two-level systems, both the 
linear decrease of the sound velocity and the linear increase of the absorption 
up to the relaxation maximum, are quantitatively accounted for by our theory. 



PACS. 61.40- amorphous and polymeric materials 
PACS. 63.50 - disordered solids, vibrational states 
PACS. 77.22G - relaxation phenomena, dielectrics 
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Low temperature properties of glasses below 1 Kelvin [[|] are satisfactorily explained 
by the assumption of localised tunnelling states (TS) with a wide distribution of energies 
and relaxation times; these TS are commonly described by a mapping on two-level systems 
(TLS) JJ. Usually only the direct (one-phonon) relaxation mechanism of TS with phonons 
is considered 0]. For this reason there is poor agreement with experiment at temperatures 
above a few Kelvin @J5|||; especially the linear increase of the absorption up to the relaxation 
peak at about 50 K and the linear decrease with temperature of the sound velocity seem 
to be a universal characteristic of amorphous substances |4],|5]||. There are attempts to 
explain that linear temperature variation by thermally activated processes ][7],[5] , or by elastic 
anharmonicity of the lattice 0, or by a modification of the standard distribution function 
for the tunnelling parameters f|. 

In this letter we provide an explanation for the temperature variation above a few Kelvin 
which relies on the two-level description and the standard distribution function, thus avoiding 
introduction of new parameters. We find in this temperature regime incoherent tunnelling 
rather than coherent oscillations; as a result, all TLS, even the symmetric ones, contribute 
to relaxation and therefore yield a more pronounced temperature dependence of sound and 
microwave propagation. In this paper we extend our previous treatment of symmetric TS 
10] to the biased case. 



The model is described by the spin-boson Hamiltonian [10] 




H = -^hA a x + ^hAa z + 7 ea z + H B , (1) 
where A denotes the tunnelling amplitude, A the bias, 7 the deformation potential, and 

(2) 

the distortion of the lattice. We consider the coupling of the TS to three-dimensional 
acoustic phonons (uj, = kv ) described by H B = J2k ^kblfik and bosonic operators fulfilling 
[6^, b\,} = 5k t k>- The model is specified by the spectral density which in Debye approximation 
is given by 

J (u) = n s (v-u k ) = 7 2 ^ 3 exp(-u/u D ) , (3) 

h 2m k v 2 

where 7 2 = 7 2 / (7i 2 gv 5 h), and u D is the Debye frequency, g the mass density and v the sound 
velocity. In the tunnelling model the parameters A and A are assumed to be distributed 
according to P(A , A) = P/A , which is equivalent to 

P 

P(e, r) de dr = -== de dr (4) 

2ry 1 — r 



with a constant P and new parameters r = A^/e 2 and e = y Aq + A 2 , where r min < r < 1. 

All dynamical information is contained in the symmetrized two-time correlation function 
C zz (t), which is calculated in the framework of the Mori-Zwanzig projection formalism [IT 
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using a mode- coupling approximation p|J12|. With the projector V = Y^ a \ a a){&a.\ =1 — Q 
and the scalar product {A\B) = ti[p eq {l/2){AB + BA)}, the equilibrium density matrix p eq = 
exp(— /3if)/tr(exp(— 0H)), the resolvent matrix C a p{z) = (er Q |[£ — z] -1 ^) {a = x,y,z) of 
the Liouvillian C* = [H, *]/H can be written as 



Here Q 



a/3 



[ z5 a/3 - Q a/3 + M aP (z) ] C M {z) 
(cr a |£|<7g) is the frequency matrix and 



-5, 



(5) 



M, 



a/3 \ 



QCa a \ [QCQ — z}' 1 \ Q£ap 



( ^yy( z 

xt 

V o 



-Eyx(z) 







(6) 



is the damping matrix with the spin-phonon resolvent Y; a p(z) = (ecr a \ [QCQ — z] 1 \eap) 
In mode-coupling approximation the memory functions are decoupled according to HJlT 

E» = CM * J(u) , 

for a — x, y, z, a, s, where we have defined the bath spectral function J {of) ■ 
coth.{phu / 2) , the weighted convolution integral 



(7) 
J{oj) 



g{u) * h{u) 



duo' cosh(/3^/2) 



2 cosh(/mc//2) cosh(/3^(cj -u')/2) 



g{uj — uj') h{ui') 



(8) 



and the resolvent functions C a {z) = C aa for a = x,y,z, and C a {z) = —i {C xy {z) — C yx {z)) 
and C a {z) = C xy (z) + C yx {z). Here the imaginary parts of the resolvent functions, i. e. the 
spectral functions, are indicated by a double prime; the real parts are obtained from these 
via a Kramers-Kronig relation. By noting C'{uf) = {uj / A ) 2 C"{uu) eq. (|5| - 0) get closed 
and can be solved numerically by iteration. They show a transition from coherent tunnelling, 
where C"{oj) has resonances at uj ~ ±e and to = 0, to incoherent tunnelling motion, where 
the three resonances have merged in one single resonance at uj = whose width narrows 
with further rising temperature. In both asymptotic regimes an analytic solution of eq. (|5| 
- 0) is possible. 

(z) In the coherent or weak-coupling regime first Born-approximation is reliable; after 
replacing in (0) C'^{uj) by the free spin-spectral function, and also discarding the spin-phonon 
interaction in p eg , one easily derives the well-known results (cf. JIB] ) 



C" z {uf) = n(a z ) 2 5{uj) , 



+ 



1 



ri 



cosh 2 (/%e/2) uj 2 + Tj 
with the usual one-phonon rate 

Ti = 2T 2 = 



r 
2 



r. 



(uJ-eY + Yl 



{uJ + ef + Yl 



71 



r'-fe* coth(J3he/2) 



(9) 



(10) 



{ii) In the incoherent or strong-coupling regime the full dynamics in C"{oj) and p eq is 
kept and the eq. (|5| - 0) are treated self-consistently. Off-diagonal correlations like T, xy {z) 



3 



and spin-polarisations (cr x ), (a z ) are now negligible. The relevant singularities of C z {z) are 
at z± = ±Q — iT 2 , Zq = —iTi where Q ps A, T 2 ~ V, Ti ~ Aq/T asymptotically. Here 
we have identified Ti x (z) with £ y (z) and have defined T by i T = T, x (z = 0). For f 3> e 
the oscillating poles have zero residue, so that the asymptotic form of the relevant spectral 
function reads 

Inserting this in the mode-coupling integral (0) yields 

For low-frequency acoustic experiments on glasses, only the relaxational pole at u = is 
relevant. For this pole the following formulae reasonably interpolate between the behaviour 
in the coherent - |TD|) and the incoherent regime ([II] - |T^) 

reiv ; r 2 + e 2 uj 2 + rf v ; 

with the relaxation rate 

I\ = J^- =: ^- , (14) 

r 2 4- e 2 t ■ 

i- r c ' min 

7min being independent of r, and 

f 7 2 e 3 coth(^e/2) =: f lph T <T 



^7 (^Y =: f M c T > T* 



(15) 



Here T* is the temperature were all thermal TLS {he < k B T) are overdamped. The condition 
Fmc(^*) = ksT*/h yields the transition temperature 

1h 1 

T* = 4f4. (16) 

Internal friction and variation of sound velocity are given by |1J 

Q^ = f-^f{u) (17) 
gv 1 

5v 1 7 2 
v 2 gv 



a X'H ; (18) 



the absorbative part }("{oj) of the dynamical susceptibility is related to the fluctuating 
part of the spectral function by the fluctuation-dissipation theorem. The bar denotes 
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the average over tunnelling systems with respect to 
e =max(kBT/h, T) one finds for the internal friction 

871-C7 2 fc| T 3 



By cutting the e-integration at 



Q- 



1 2L /-» ?tf M c ( E 
~ 2 ^ 2fc s T V2 



arctan 



^ r M c r min ^ ~j 



T > T* 



(19) 



and - by adding the contribution of the resonant part Sv/v\ 
of the sound velocity 



Cln(T/T ) - for the change 



Sv 



-IClnJ 



T <T* 



T 1 * << T 1 << T 1 * 

res ^ ^ 



(20) 



16 



ln( 



r 2 



T*huj/k L 



with C = P^/gv 2 . T* es separates the regimes where the resonant (T < T* cs ) and the 
relaxational process (T > T r * s ) prevails. Below T* one finds the well-known logarithmic 
temperature dependence of the sound velocity and the constant internal friction. At T = T* 
the temperature dependence changes to a linear increase in the absorption and a linear 
decrease in the sound velocity. 

From a recent experiment on Suprasil W [14j] one finds a transition from the plateau to 
the linear increase at about T* w 6 K which corresponds according to to a deformation 
potential of about 7^2 eV. Here we have used C = 2.8 x 10~ 4 [14], g = 2.2 g/cm 3 , 
vi = 5.8 x 10 5 cm/sec, v t = 3.75 x 10 5 cm/sec and t| w 27 t 2 with 7 2 /f 5 = '-)j/vj + 2^ 2 /v 5 . 
With these values we calculate the slope with respect to temperature; for a comparison with 
experiment see Table 1. Both the prefactor and the logarithmic variation of the sound 
velocity with frequency show full agreement. 

In Fig. 1-2 we have plotted our theoretical results together with the experimental data for 
Suprasil W |14|. At temperatures between 100 mK and 15 K there is full agreement between 
experiment and theory; for both absorption and sound velocity the measured data could be 
reproduced with the same numerical values for 7 and C. A similarly good agreement has 
been found at other frequencies and in recent experiments on GeC>2 [[15[ . 

Finally, we comment on the absorption peak at about 30 K. According to (19) we find in 
our theory T max oc y/uj. However, experimentally the frequency dependence of the relaxation 
maximum is found to be much weaker 



nr 



This would indicate the onset of thermally 
activated processes to occur at some temperature below T max which yields a logarithmical 
variation with frequency, rather than the square root dependence. For that case the quoted 
value of r min in Fig. 1 has no physical relevance. We stress that this does not affect the 
temperature variation below T max . 

In summary, we have shown that the experimental data in the absorption and the sound 
velocity up to the relaxation peak can be explained in the framework of the tunnelling 
model. The novel features arise from the incoherent dynamics of the tunnelling motion 
at temperatures above T*. In particular, this provides a new relaxation mechanism which 
accounts well for the experimental findings in glasses above 5 K. 
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CAPTIONS 



Fig. 

HI. 



1: Internal friction for Suprasil W at 11.4 kHz. The data (o) are from Classen et al. 
In our theory (-) we have used C = 2.8 x 1CT 4 and 7^ = 2.2 eV 

2: Relative variation of the sound velocity for Suprasil W at 11.4 kHz. The data (o) 

we have used the same numerical value for C 



Fig. 

are from Classen et al. 
and 7^ as in Fig. 1 



14 . In our theory 
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TABLES 



Table 1: Comparison of theoretical and experimental results for Suprasil W in the incoherent 
regime T > T*. The experimental values are taken from [|14 





Experiment 


Theory 


Q- 1 


(7 ±2) x l(r 5 T/K 


6.1 x 10" 5 T/K 


Sv/v 


2 x 10- 5 T/K m(1.6 x 10- 13 cj/2vrHz) 


2.2 x 10- 5 T/K \n(T*htu/(k B T 2 )) 
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